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Abstract. The Landau Hamiltonian governing the behavior of 
a quantum particle in dimension 2 in a constant magnetic field is 
perturbed by a compactly supported magnetic field and a similar 
electric field. We describe how the spectral subspaces change and 
how the Landau levels split under this perturbation. 
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1. Introduction 

The paper is devoted to the study of the spectrum of the Schrodinger 
and Pauli operators in the plane, with nonzero constant magnetic field 
perturbed by smooth compactly supported magnetic and electric fields. 

The Landau Hamiltonian describing the motion of the quantum par- 
ticle in two dimensions under the influence of the constant magnetic 
field is one of the classical models in quantum physics. It was intro- 
duced and studied in late 1920-s, see [10]. It turned out that the system 
possesses rather unusual spectral properties: the spectrum consists of 
discrete eigenvalues lying at the points of an arithmetic progression and 
all these eigenvalues have infinite multiplicity. These eigenvalues are 
traditionally called Landau levels (LL) and the corresponding spectral 

subspaces are called Landau subspaces. 

l 
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A natural question arises, what happens with the spectrum of the 
Landau Hamiltonian under perturbations. If the perturbation is strong, 
say, it is a Coulomb electric field, it turns out more reasonable to con- 
sider the constant magnetic field as a perturbation, and it leads to 
Zeeman effect. On the other hand, one can consider weak perturba- 
tions of the Landau Hamiltonian by a fast decaying or even compactly 
supported electric or/and magnetic field. In this case it is natural to 
expect that the Landau levels will split, and the question consists in de- 
scribing this splitting as well as the behavior of the spectral subspaces 
under the perturbation. 

The case of the perturbation by an electric potential, moderately, 
power-like decaying at infinity, was first studied by Raikov in [12]. It 
was shown that each LL splits into an infinite sequence of eigenvalues, 
converging to the originating LL. The rate of convergence of these 
eigenvalues is also power-like and it is determined by the rate of decay 
of the perturbation at infinity, the asymptotics being described by a 
semi-classical formula. Advanced methods of micro-local analysis were 
used in [12]. 

The case of a fast decaying (or compactly supported) electric po- 
tential was dealt with much later, in [13], see also [11]. It turned 
out that the Landau levels still split, however in an unusual, not a 
semi-classical way. Most profoundly this can be seen for a compactly 
supported perturbation. The semi-classical considerations, formally 
applied, would lead to a finite number of eigenvalues splitting away 
from the LL. However, in fact, not a finite but an infinite number of 
eigenvalues split away, but they converge to the LL extraordinarily 
fast, super-exponentially. Moreover, the leading term in the asymp- 
totics does not depend on the perturbation at all, and the next term 
(as it was found recently in [5]) is determined by the geometry of the 
support of the perturbation. This result can be also understood as 
justifying the first order Raleigh-Schrodinger perturbative correction, 
with properly chosen initial basis in the Landau subspace. Methods of 
micro-local analysis do not work here, and the results are obtained by 
a combinations of methods of complex analysis and variational consid- 
erations. 

Considerably less is known about the perturbations of the Landau 
Hamiltonian by a weak magnetic field. The complication here has 
its roots in the fact that it is not the magnetic field itself but its 
potential that enters in the quantum Hamiltonian. Unless the total 
flux of the perturbing magnetic field vanishes, the perturbation of the 
(Schrodinger or Pauli) operator contains rather slowly decaying coeffi- 
cients, so it is fairly strong even for a compactly supported perturbation 
of the field and it even may be not relatively compact if the perturba- 
tion goes to zero at infinity not sufficiently fast. Iwatsuka [9] proved 
that invariance of the essential spectrum still takes place, so Landau 
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levels are the only possible limit points of eigenvalues in the gaps be- 
tween them. After that, there were very few results in the spectral 
analysis of this problem (apart from the rather specific situation with 
the lowest LL). In the paper [4] the perturbation by a single Aharonov- 
Bohm solenoid was considered. For the standard self-adjoint extension, 
by means of explicitly solving the eigenvalue problem, it was established 
that the Landau levels are still intact, the eigenspaces change in a con- 
trolled way, and a finite number of eigenvalues split away from each 
LL. Hempel and Levendorski in [8] and Besch in [3] investigated the 
behavior of the eigenvalues in the gaps when a large parameter tending 
to infinity stands in front of the perturbation. The number of eigen- 
values between the Landau levels grows quasi-classically, as a power 
of the parameter, but nothing was found out about the distribution 
of the eigenvalues for a fixed value of the parameter. However it was 
shown that one should not expect monotonicity of the eigenvalues in 
this parameter, even if the perturbing field has constant sign. 

The aim of the present paper is to find out what happens with the 
spectrum and spectral subspaces of the Landau Hamiltonian under 
smooth compactly supported perturbations of the magnetic field, with 
additional perturbation by a compactly supported electric field (thus 
both the case of the magnetic Schrodinger and of the Pauli operators 
are taken care of). We show that the Landau levels split superexpo- 
nentially, and under certain positivity conditions find the asymptotics 
of eigenvalues as they approach the LL. We find also a rather exact 
approximation to the spectral subspaces of the perturbed operator. 
It turns out that the spectral subspaces are perturbed fairly strongly. 
This effect is absent for the perturbation by the electric field only: we 
show that in the latter case the spectral subspaces change very weakly, 
in proper terms (the behavior of spectral subspaces was not investi- 
gated in [13], [11]). 

Our approach is based upon the study of approximate creation and 
annihilation operators. The commutation relations for these operators 
(the deformed Heisenberg relations) enable us to construct approximate 
spectral subspaces. It is known that for the unperturbed operator the 
Landau subspaces are obtained by repeatedly applying the creation 
operator to the explicitly described infinite-dimensional space of zero 
modes of the Pauli operator. The Pauli operator with the perturbed 
magnetic field possesses an infinite-dimensional space of zero modes as 
well, moreover, this subspace can also be explicitly described. It turns 
out that the application of the same procedure for the perturbed oper- 
ator, but using the approximate creation operator instead, produces a 
very good approximation for the spectral subspaces for the perturbed 
operator. 

In the papers [13], [11], [5], the splitting of Landau levels under 
the electric perturbation was associated with the spectrum of certain 
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Toeplitz-type operators acting in the spectral subspaces. In our case, 
the study of the spectrum of the perturbed operator is based upon the 
analysis of Toeplitz-type operators in the space of zero modes, general- 
ized in the sense that a differential operator replaces the multiplication 
by a function in the Toeplitz expression. 

We prove the asymptotics of the eigenvalues split from the Landau 
level under a certain positivity condition imposed on the effective po- 
tential reflecting the combined influence of the electric and magnetic 
perturbations upon the spectrum. If this condition is not fulfilled, we 
can only prove the superexponential estimate from above. The obstacle 
here is the lack of knowledge of the spectrum of the Toeplitz-like op- 
erators with indefinite weight in the space of zero modes, which turns 
out to be a separate hard problem in complex analysis. 

In Section 2 we define our main operators and establish a series 
of relations between them. Further, in Section 3, we construct the 
approximate spectral subspaces. In Section 4 we prove that they form 
a very good approximation for the actual spectral subspaces. As a by- 
product, we show that the spectral subspaces 'almost' do not change 
if just an electric perturbation is present. Section 5 is devoted to the 
analysis of the Toeplitz-like operators. In Sect. 6 we establish the 
estimates, and when possible, the asymptotics in the Landau levels 
splitting. The last section contains the proofs of several rather technical 
formulas and estimates. 

Different results involve different smoothness conditions imposed on 
the perturbations. The upper estimate for the eigenvalue splitting is 
proved for C 2 compactly supported perturbations. On the other hand, 
the lower estimates granting the eigenvalue asymptotics, as well as the 
result on the approximate spectral subspaces are proved for infinitely 
differentiable perturbations. These latter smoothness conditions can 
be replaced by some finite smoothness, depending on the Landau level 
under consideration, but we do not go into these details. 

The conditions of the perturbations to have compact support can 
be replaced by the requirements of controlled decay at infinity. In 
that case the results on the eigenvalue splitting and on the spectral 
subspaces approximation can also be obtained, with the eigenvalue 
estimates and asymptotics as well as the quality of spectral subspaces 
approximation depending on the rate of decay. These resulst will be 
published elsewhere. 

The second author was partly supported by a grant from the Royal 
Swedish Academy of Sciences, he also thanks Chalmers University of 
Technology in Gothenburg for hospitality. 

2. Magnetic Schrodinger and Pauli operators 

We will denote the points in the plane R 2 by x = (x 1 ,x 2 ); it is 
convenient to identify M 2 with C by setting z = xi+ix 2 . The Lebesgue 
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measure will be denoted by dx and the derivatives by dk = d Xk ; we set, 
as usual, 3 — (<9i + id 2 )/2, d — (di — id 2 )/2. 

The constant magnetic field will be denoted by B°; we fix the sign, 
B° > 0. With this magnetic field we associate the magnetic potential 
A°(x) = (Al,A$) = ^-(—X2,xi). Then the (unperturbed) magnetic 
Schrodinger operator in L 2 (IR 2 ) is defined as 

S)° = -(V -iA°) 2 . (2.1) 

The Pauli operator describing the motion of a spin-| particle acts in the 
space L 2 (1R 2 ) 2 of two-component vector functions, it has the diagonal 
form, = diag <P° ), where = #° ± B°. 

The spectrum of the Schrodinger operator was found by Landau in 
1928 (see [10]) by separation of variables. However more conveniently 
the spectral structure of these operators is described by means of cre- 
ation and annihilation operators introduced by Fock [6] and reflecting 
the algebraical structure of the problem (see, for example the detailed 
exposition in [7] or [2]). Let us introduce the complex magnetic poten- 
tial A° = A\ + iA° 2 and set 

£3° = -2id - A°, 5° = -2id - A°. (2.2) 

These operators can be also expressed by means of the scalar potential 
of the magnetic field, the function 

= ^-\z\ 2 (2.3) 

solving the equation = B°: 

= -2ie^de^\ ;Q° = -2ie-*°de*°. 

The creation and annihilation operators £3°, £3° satisfy the following 
basic relation 

[0°,£F] =2B°. (2.4) 

This relation, meaning that the operators H°,0°, 2B° generate a rep- 
resentation of the Heisenberg algebra, together with the formulas 

enable one to describe the spectrum completely, and the construction 
is presented in numerous physical and mathematical sources. Since 
0° = £3°*, the equation u = is equivalent to 

;Q° M = e-*°d{e*°u) = 0. 

This means that / = e r M is an entire analytical function in C, such 
that after multiplication by e~*° it belongs to L 2 . The space of such 
functions / is called Fock or Segal-Bargmann space T (see [1], [2], 
[7] for extensive discussion); when we need to specify the particular 
weight involved in the definition of T we will use the notation JF B o. 
So, the null subspace of the operator i.e., its spectral subspace 
corresponding to the eigenvalue A = 0, is C = e-~^° T ■ After this, 
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the commutation relations (2.4), (2.5) imply that Ai = 2B° is the 
lowest eigenvalue of the operator %$° + = + 2B°, with the same 

spectral subspace, and thus C\ = 0°£o is the spectral subspace of 
*J}° with eigenvalue Ai and so on. So, the spectrum of consists of 
eigenvalues, Landau levels A q = 2gB°, q — 0, 1, ... , and the spectra of 
consist, respectively, of A q + B° and A q + 2B°. The operators 
0°, 0° act between Landau subspaces C q = £l oq Co, q = 0, 1, . . . , 

& : C q h+ C q+1 , Q°:C q ^ C q -!, 0° : C h-> 0, (2.6) 

and are, up to constant factors, isometries of Landau subspaces. 

It is also convenient to introduce the spectral projection P° of ^3° 
corresponding to the eigenvalue A q = 2gB°. The above reasoning can 
be expressed as 

P° = C q - 1 &' 1 P°£l oq , C q = q\{2B ) q . (2.7) 

Now we introduce the perturbation b e C^°(1R 2 ) of the magnetic 
field, and we set B = B° + b. Let ip be a scalar potential for the 
field b, a solution of the equation = b. Of course, tp is defined up 
to a harmonic summand, the choice of corresponds to the choice of 
gauge. The magnetic potential a = (ai,a 2 ) = (— <9 2 t/>, <9i ■?/>), curia = b, 
is thus determined up to a gradient, and the complete scalar and vector 
magnetic potentials are 

$ = r + ^, A = A° + a. (2.8) 

The crucial observation here is that although b has compact support, 
we should not expect the same, or even a rapid decay, from the mag- 
netic potential a, unless the total flux of b, J b(x)dx is zero. This can 
be seen in the easiest way by applying Green's formula, saying that the 
circulation of the vector field a along the circle with radius R equals 
the flux of b = curia through the disk \x\ < R. Thus, if the total flux 
of b is nonzero and R is taken large enough, so that the support of b 
lies in the disk, the circulation of a is bounded away from zero, and 
thus a may not decay faster than CR _1 , at least in some directions. 
This corresponds to a logarithmical growth rate of ip. 

We define the perturbed magnetic Schrodinger operator as in (2.1), 
with A° replaced by A: 

£ = _(V- ? A) 2 , 
and the components of Pauli operator as 

<P± = Sj ± B. 

It is easy to observe that the difference between ^ and ,fj contains an 
operator of multiplication by A° • a, the latter function does not decay 
at infinity (again unless the total flux of b vanishes), and thus looks like 
being not a relatively compact perturbation of S)°. However, thanks to 
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the special form of this term, the perturbation is still relatively compact 
if a — > at infinity, as was noticed by Besch [3] . 

Now let us look at the algebraic structure related with the perturbed 
operators. The perturbed creation and annihilation operators are de- 
fined similarly to (2.2): 

= -2id - A, = -lid - A, 

where A is the complex magnetic potential, A = A° + (ai + ia 2 ). The 
representation via the scalar potential takes the form 

= -2ie*9e-*, = -2i e -*de*. 

Most important for us are the commutation relations for 0, 0: 

[0,0] = 2B = 2B° + 2b, (2.9) 

together with 

<P+ = 00, «p_ = 00, - = 2B = 2B° + 2b, (2.10) 

£ = 00-B = 00 + B. (2.11) 

The relations (2.9)-(2.11) look similar to (2.4), (2.5), however the 
crucial difference is that the new ones contain variable functions on 
the right-hand side, do not represent the Heisenberg algebra any more 
and thus cannot be used for the exact spectral analysis. The spectra of 
Schrodinger and Pauli operators do not determine each other any more. 
The only information that one can obtain immediately, is the descrip- 
tion of the lowest point of the spectrum of Since, again, = 0*, 
the equation = is equivalent to 0-u = 0, or d(exp($f)u) = 0. 
So the function f = u exp \& is an entire analytical function such that 
u = exp(— ty)f G L 2 . The space of entire functions with this property 
is, obviously, infinite-dimensional, it contains at least all polynomials 
of z, although it does not necessarily coincide with Fock space, since 
exp(— ip) may have power growth (or decay). We denote this space by 
jF b . 

The same reasoning, applied to shows that the null subspace 
of 0, as well as the one of is zero: it consists of anti-analytical 
functions which belong to L 2 after the multiplication by the growing 
exponent exp(\&), and there are no such functions, except zero. 

We denote the null-space of the space of zero modes, by H,q. In 
[14] rather general sufficient conditions were found for infiniteness of 
dimension of Tto', in particular they even do not require the magnetic 
field to be a weak perturbation of the constant one. 

As it follows from the relative compactness of the perturbation, by 
Weyl's theorem, the essential spectrum of the perturbed operator 
consists of the same Landau levels A q , and the eigenvalues in the gaps 
and below A may only have A q as their limit points. This latter fact 
was established much earlier by Iwatsuka [9], and the reasoning based 
upon the commutation relations remains valid even in the case of a 
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slowly decaying magnetic perturbation, when the perturbation of the 
operator is not relatively compact any more. Since our considerations 
in the paper can be considered as an extension of the ones of Iwatsuka, 
we reproduce the reasoning from [9]. 

Theorem 2.1. ([9]) The essential spectrum of the operator consists 

of the points A q , q = 0, 1, The essential spectra of $) and 

consist, respectively, of the points A q + B° and A q + 2B° . 

Proof. Note first that the operator of multiplication by b is relatively 
compact with respect to S) and ?fi±. Thus, by Weyl theorem, the oper- 
ators and + 2B° = — 2b have the same essential spectrum, 
i.e. the essential spectrum of ^3+ is just the one of shifted by 2B°. 
So, ^3+ has no essential spectrum in [0, 2B°). On the other hand, the 
nonzero spectra of and coincide, since these operators are just 
products of and in the opposite order. Therefore has no essen- 
tial spectrum in (0,2B°), and 2B° belongs to the essential spectrum. 
Again, this means that the interval (2B°, 4B°) is free from the essential 
spectrum of and 4B° belongs to it. By repeating this reasoning, 
the complete description of the essential spectrum is obtained. □ 

Now we add a perturbation by the electric potential. Let V{x) be a 
real valued function in C^°(M 2 ). We introduce the operators 

S)(y) =S) + v, qj±(y) = <p± + v. (2.12) 

Since the operator of multiplication by V is relatively compact with 
respect to 5},^P±, the operators (2.12) have the same essential spectra 
as the respective unperturbed ones. In this paper we are going to study 
the behavior of the eigenvalues of S)(V),^P±(V) near these Landau 
levels. 

3. Approximate spectral subspaces 

In this section we construct the approximate spectral subspaces of 
the operators (2.12) under the condition b e C^°. To do this, we need 
some more properties of the space Ho of zero modes of ^P_. We must 
keep in mind that the ambiguity in associating the vector potential 
with the magnetic field leads to the ambiguity in defining the magnetic 
operators themselves. The gauge transformations of the form u i— > e te u 
with real functions Q(z) such that e l@ is smooth, realize the unitary 
equivalence of magnetic operators constructed by different choice of the 
potential. 

As it was explained in Section 2, the space Ho consists of functions 
u(z) having the form u(z) = e _ *^/(z) where f(z) is an analytical 
entire function such that e~^ f G L 2 , and ^(z) is the potential in 
(2.8). Note first that the choice of * does not affect H , up to a gauge 
transformation. In fact, let $ be a harmonic function in R 2 , so that 
ty' — + $ is also a potential for the magnetic field B = B° + b, 
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A\l/' = B. Denote by G the function harmonic conjugate to $, so 
that F = $ + i<d is an entire analytical function. This gives us the 
representation 

u(z) = e- ie ^e-*'( z \e F(z) f(z)), 

so the multiplication by e 1 ®^ is the gauge transformation in question. 

However the perturbation of the magnetic field leads to the pertur- 
bation of the space of zero modes that cannot be reduced to a gauge 
transformation. In fact, if there were a real function Q(z) such that 
e l@ Ti. = Cq, this would mean that the function e~^ +l0 is analytical. 
However in this case the function —ip, as the logarithm of the absolute 
value of an analytical function, must be harmonic, which it is not. 

Nevertheless there is a close relation between Ho and C , not re- 
duced to a gauge transform. We fix the potential ip, defining it as the 
logarithmic potential of the field b, 

rl>(x) = (27T)- 1 J b(y) log \x - y\dy. (3.1) 

The behavior of the potential ip at infinity is determined by a, the 
total flux of the perturbing field divided by 2n, 

a = (27T)- 1 J b{x)dx. 

If a — 0, is a bounded function, as well as e^. So, the multiplica- 
tion by the bounded function m = maps Co to Ho- This mapping 
is, obviously, invertible. 

Now suppose that o > 0. Then the function ip(x) grows at infinity 
as cr log | re j , m(x) = grows as and it is not bounded any 

more. The function e~^^ is bounded. Thus the multiplication by 
mr 1 is a bounded operator from Ho to £ - The inverse mapping is 
not bounded, but we show that, nevertheless, it is well defined on a 
certain subspace in Co- Denote by v the smallest integer larger than 
a. Fix v arbitrary different points z±, . . . ,z v in C. Denote by p(z) the 
monic polynomial with simple zeros at the points z±, . . . , z v and by Co 
the space of functions in C having zeros at these points. Since the 
evaluation at a point is a bounded functional in the Fock space, Cq is 
a closed subspace in Cq, with codimension v. Now consider the space 
Cq = p(z)' 1 Co. By our construction, this space is still a part of Co 
(however is not a closed subspace any more). Since at infinity we have 
e 4>(z) _ 0(p(z)~ 1 ), the operator of multiplication by maps Cq to 
Ho- In the case of a negative a, the spaces Tto and Co change places in 
the above reasoning. 

The space 7i is closely related to the spaces $ p introduced in [1], 
Section 3, consisting of entire analytical functions f(z) that belong to 
L 2 after being multiplied by 9 p (z) = (1 + \z\ 2 )~ p / 2 exp(-(B°/4)|z| 2 ). 
Due to (3.1), the function £(z) = ip(z) — f log(l + \z\ 2 ) is bounded, so 
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we can write 

e -*W = 0_ ff (z) e -«*>. 

Therefore our space 7i can be obtained from Bargmann's space Q_ a (z)$ a 
by multiplication by the bounded and boundedly invertible function 
e -?0). This means that many properties of this Bargmann space are 
inherited by our space TCq. What is important for us at the moment 
is the fact established in [1], Section 3.6, that the functions of the 
form 9 cr (z)z n form a basis in 8 a (z)$ cr . Therefore this space possesses 
a dense set consisting of rapidly decaying functions (having the form 
0- er (z)p(z), with polynomial p), and thus our null subspace TCo pos- 
sesses a dense subspace of rapidly decaying functions as well. 

Let 5 q = (A q — 7, A q + 7), q = 0, 1, 2, ... , be intervals of the same 
size on the real axis, centered at the Landau levels A q = 2gB°, and not 
containing the neighboring Landau levels A q ±i. We choose the size of 
5 g in such way that neither of these intervals has the eigenvalues of 
at its endpoints. Moreover, since the lowest LL A = is an isolated 
point of the spectrum of ^P_, we can choose the size of the intervals in 
such way that So contains only this point of spectrum. We denote by 
7i q the spectral subspace of corresponding to the interval 5 q and 
by P q the corresponding spectral projection. Since, by Theorem 2.1, 
the spectrum of is discrete between Landau levels, the change of 5 q 
leads only to a finite-rank perturbation of P q . As usual, the spectral 
projection P q can be be expressed by means of the integration of the 
resolvent of along a closed contour T g in the complex plane, not 
passing through the eigenvalues of and containing inside only those 
eigenvalues that lie in 5 q . Again, using the discreteness of the spectrum 
of between the Landau levels, we can choose these contours in such 
way that all they are obtained from Tq by the shift along the real axis 
in the complex spectral plane, T q = T + 2gB°. 

Now we are going to establish several properties of the subspaces C q , 
projections P q and some related operators. 

First, note the simple fact following directly from the spectral theo- 
rem. 

Proposition 3.1. For any q = 0, 1, . . . , and any polynomial p(X) the 
operator p{ty J) P q is bounded, moreover the operator {pity ) — p{A q ))P q 
is compact. 

In fact, by the spectral theorem the nonzero spectrum of the operator 
p(^3_)P 9 consists of the points p(Xj) where Xj are all points of spectrum 
of in S q and thus all p(Xj) live in a bounded interval. Moreover, 
p(Xj) may only have p(A q ) as their limit point, by Theorem 2.1. 

The following lemma will enable us later to establish a much stronger 
compactness property. 

Lemma 3.2. Let each ofTj, j — 1, . . . , N be one of operators or 0. 
Then for some constants C, C , for any u in the domain of the operator 
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\\T{r 2 . . . T N u\\ 2 < C{<$%, u) + C"|M| 2 . (3.2) 

The proof of Lemma 3.2 can be found in Section 7. 
Now we can establish the compactness property. 

Proposition 3.3. Let T ± , . . . ,T N be a collection of operators, each be- 
ing 13 or 13 ; and let hj, j — 0, . . . N be functions with all derivatives 
bounded, T = h Tihi . . . T N h N . Then for any q and for any polynomial 
p(X) the operators T(p(*p_) — p(A q ))P q T are compact. 

Proof. By commuting functions hj and operators Tj (moving all func- 
tions to the left), we can transform the left operator T to the sum of 
terms having the form h K T K where h K are bounded functions and T K is 
a product of no more than N operators 13, 13. Similarly in the copy of 
T that stands to the right from P q , we move all functions to the utmost 
right positions, to get the representation of T as a sum of terms T^h^, 
being the product of a bounded smooth function h^ and no more than 
iV operators £3, 13. 

For each of the terms arising in this way in the decomposition of 
T(p(«p_) -p(A q ))P q T, we can write 

-p(K))P q TJix = [h K %(^ + l)- x ]x (3.3) 

[cp? + i)(ppp_) - P (A q ))p q (v N + i)]w N + irxu 

In (3.3), the first factor in brackets is bounded by Lemma 3.2, and the 
middle factor is compact by Proposition 3.1. The last factor in brackets 
is also bounded, by Lemma 3.2 applied to the adjoint operator. □ 

The result we have just proved indicates that, probably, the eigen- 
values of the operator lying near A q , must converge to A q rather 
rapidly. This will be, in fact, established later on. 

Now we are going to describe the main construction of the paper, 
the approximate spectral subspaces of the perturbed operator. Note, 
first of all, that it is sufficient to consider the operator ^3_(V). In fact, 
due to the commutation relations, we have 

Sj(V) = <p_(V + b) + B°, % + (V) = + 2b) + 2B°, 

and thus the operators differ from by a shift and by the electric 
type perturbations V + b, V + 2b. We start by finding approximate 
spectral subspaces of Adding V will be an easier task after this. 
The subspaces approximating H q will be defined in the following 

way: 

Go = Ho, Q q = 5 9 o , q = 1, 2, ■ ■ ■ ■ (3.4) 
So we mimic the construction of the eigenspaces of the unperturbed 
Landau Hamiltonian, see (2.6). Our approximate spectral subspaces 
are obtained by applying the creation operators to the space of zero 
modes. 
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Of course, since we apply the unbounded operator 0, we must show 
that all the time we do not leave the space L 2 , and moreover, we have to 
prove that the subspaces Q q are closed. Both these properties, as well 
as some other results of the paper will be based upon the important 
Proposition 3.4 to follow. But first we need a definition. 

Many times in the paper we will encounter operators that are prod- 
ucts of several copies of 0, 0, functions B, b, U (for some smooth func- 
tion U) and their derivatives. Let us assign the weight 2 to B°, B, b, U, 
the weight I + 2 to the order / derivatives of b, U and the weight 1 to 
and to 0. For a product of operators and functions the weights of 
factors are added. 

Proposition 3.4. Let q > 0. 

(1) There exists a function Z q [b] G depending only on q, B°, 
and b such that for any u G Ho, 

\\Qi u \\ 2 = C q \\u\\ 2 + (Z q [b]u,u), C q = q\(2B°) q . (3.5) 

The function Z q [b] is a polynomial of b and its derivatives up 
to the order 2q — 2 with coefficients depending on B°. All sum- 
mands in Z q [b] have weight 2q. The term with the highest de- 
rivative equals 2A q b; the term linear in b and not containing 
derivatives equals C'B o9_1 b ; C = 2 q q\q. 

(2) There exists a function X q [b, U] G depending only on q, B°, 
b ; and U such that for any u G Hq. 

(C/0%,0%) = (X q [b,U]u,u). (3.6) 

The function X q [b, U] is expressible as an order 2q linear differ- 
ential operator acting on U, with coefficients depending polyno- 
mially on b, its derivatives, and B°. All the terms in X q \b,U\ 
have weight 2q + 2. The term with the highest derivative ofU is 
A q U. The term with zero order derivative ofU is C q U +Z q \b]U . 

(3) There exists a function y q \b, U] G depending only on q, B°, 
b ; and U such that for any u G 7io- 

(UQ q u,~Q q+1 u) = (y q \b, U]u,u). (3.7) 

The function y q [b,U] is expressible as an order 2q + 1 linear 
differential operator acting on U , with coefficients depending 
polynomially on b, its derivatives, and B°. All the terms in 
y q \b, U] have weight 2g+3. The term with the highest derivative 
of U is (—2i)A q dU. The term with zero order derivative of U 
is C q U + Z q [b]U. 

Proof. First, in order to make the main mechanism of the proof more 
clear, we derive the relation (3.5) for q — 1 and q = 2. It is sufficient 
to prove (3.5) for u = exp(— ^)p(z) where p(z) is a polynomial: such 
functions are dense in Ho, and then we can pass to general u by the 
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closure. In the process of our, rather lengthy, calculations we will 
repeatedly use the commutation relations (2.9) as well as 

[0, h] = -2idh, [0, h] = -2idh, for h e C°°. (3.8) 

So, let us consider the expression ||0-u|| 2 , u e Tio- We have 

||0w|| 2 = (5u,5u) = (QQu,u) (3.9) 

= (2B° + 2b + 00w, w) = 2B°\\u\\ 2 + (2bu,u) 

(on the last step we used that 0w = 0). Thus we have Zi[b] = 2b. 
The case q = 2 requires some more calculations. We start with 

\\Wu\\ 2 = (£l 2 Q 2 u,u). 

Then we are going to commute and so that only 0, and not 
would act on u, so we will move all copies of to the left and all copies 
of to the right: 

(0 2 0V u) = (0(2B° + 2b + 00)0m, u) = 2B°||0«|| 2 (3.10) 

+2(0b0w,w) + (0000w,w). 

The first term in (3.10) has already the form we aim for, by (3.9). In 
the second term we commute and b using (3.8): 

2(0b0w, u) = 2(b00w, u) + 2(-2i(9b)5«, u). (3.11) 

In the first term in (3.11) we commute and 0, and in the second 
one commute — 2i(dh) and 0: 

2(0b0w, u) = 4(b(B° + b)u, u)+ 

2(b00w, u) + 2(S(-2i9b)u, u) + (8(ddb)u, u). 

Since Q.u = 0, we arrive at 

2(0b0 M , u) = 4(b(B° + b)u, u) + 2((Ab)u, u). 

In the last term in (3.10) we commute the first pair and and then 
the second one: 

(0000m, u) = (00 2 0m, u) + 2B°(00m, u) + 2(b00«, u) (3.12) 

= 2B°(0m, 0m) + 4(b(B° + b)u, u) + 2(b00w, u). 

On the right-hand side of (3.12), the first term is again calculated by 
(3.9), the second term has the required form, and the last one vanishes. 
Now we collect all the transformed terms in (3.10) and obtain 

||0 2 w|| 2 = 2B°(2B°|H| 2 + (2bu, «)) + 4(b(B° + b)u, u)+ 
2(Abw, u) + 2B°(2B°||m|| 2 + (2bu, u)) + 4(b(B° + b)u, u) 
= 2(2B ) 2 \\u\\ 2 + (Z 2 [b]u,u). 

with 

Z 2 [b] = 8b 2 + 12B°b + 2Ab. 



14 



ROZENBLUM AND TASHCHIYAN 



Now we can formally describe how the procedure goes in the general 
case. In the expression 9 <? we start commuting and 0, as well as 
these operators with the functions arising by the previous commutings. 
All the time we keep moving to the right and to the left. As soon 
as reaches the utmost right position or reaches the utmost left 
position, the corresponding term vanishes. We do this until no oper- 
ators 0,0 are left, and only a function C q + Z q [b] remains. To find 
C q and describe the structure of £ ? [b], we need only to do some book- 
keeping of the terms arising when the operators and/or functions are 
commuted. The formulas (2.9), (3.8) show that the weight is preserved 
under commutation. So the constant term C q has to be proportional to 
B 09 , and it can arise only when the operators 0, are commuted, not 
the operators with functions. The number of such commutings needed 
to re- arrange all terms equals q\. The term linear in b and not con- 
taining derivatives results from the same procedure and has the form 
CgB° 9_1 b, C' q = 2 n q\q. The function Z q [h] arising as a result of moving 

all 0-s to the right and all 0-s to the left will be a polynomial of b 
and its derivatives with weight 2q. The highest derivative of b present 
in Z q [h] will thus be of order 2q — 2, and this term will be linear in b. 
Such term arises when 2B appears after the first commutation of 0, 



and then q — 1 copies of pass 2b on their way to the right and q — 1 
copies of pass it on their way to the left, each such pair producing 
the Laplacian acting on b. 

The above reasoning takes care of (3.5). Parts 2) and 3) of the 
proposition are proved in the same way, with a similar bookkeeping in 
order to determine the structure of the resulting functions. We just 
demonstrate how Part 2) works for q = 1: 



(QUQu,u) = (QWu,u) + (0[U,Q]u,u) = (QWu,u)(= 0)+ 
(2BUu,u) + ([l/,5]0u,u)(= 0) + ([0, [U,~5\u,u) = (^[b, U]u, u); 



When applying Proposition 3.4 and similar results, we need a certain 
compactness property. Such facts were used persistently in [13], [11], 
[5], but for the case of a constant magnetic field only. 

Lemma 3.5. Let W(x) be a bounded function with compact support. 
Then for any differential operator £ with bounded coefficients the qua- 
dratic form 



in 0^09, 



0"0 9 = 0^ 1 000 9 + 9 ~ 1 (2B° + 2b)0' 



,9-1 



X 1 [b,U] = 2BU + AU. 



□ 




is compact in the space Hq. 
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Remark 3.6. In section 5 we will find much more detailed information 
about such forms. However the proof we give here is also useful. We 
return to this approach later. 

Proof. Let u = e~^ / be a function in Tl , so f(z) is an analytical 
function. 

Take some R such that the support of the function W lies inside the 
circle Cr, with radius R centered in the origin. For each r e (R, 2R) 
we write the Cauchy representation for /: 

f(z) = (2m)- 1 f fiOiz-CY'dC (3.13) 

Then for some fixed function £(r) G C °°(R,2R), / £{r)dr = 1, we 
multiply (3.13) by £(r) and integrate in r from R to 2R. This gives 
the integral representation of f(x), \x\ < R, in the form f(x) = 
f R< ^ <2Ii K(x,y)f(y)dy, with smooth bounded kernel K(x,y). After 
applying £ in x variable, we obtain a similar integral representation for 

£u(x)= [ e- nx) K st (x, y y { y ) u(y)dy=(& z u){x). 

jR<\y\<2H 

The integral operator IVF) 1 / 2 .^ has a bounded kernel with compact 
support and therefore is compact in L 2 , and thus the quadratic form 
to[u] can be written as tt)[u] = (sign WlWl 1 / 2 .^^ |H^| 1//2 ^ £ -u) and 
therefore it is compact. □ 

Now we are able to justify our construction of the spaces Q q . 

Proposition 3.7. The sets Q q defined in (3.4) are closed subspaces in 

Proof. The fact that Q q C L 2 follows directly from Proposition 3.4. 
Next, the relation (3.5) can be written as 

(P o 9 5%, u) = C q {u, u) + {P Z q [b]u, u); u = V ue H (3.14) 

(recall that Pq is the orthogonal projection onto Ho, PqU = u). In the 
second term in (3.14), by Lemma 3.5, the operator P ^g[b] is compact 
in Ho, and therefore we can understand (3.14) as showing that the 
operator C^PoH 9 is a left parametrix for £l q : H — > L 2 . This implies 
that the range of £l g is closed. □ 

The null space of £3 and therefore of £l q is zero. Consider the operator 
£l g as acting from Q = Tio to Q q . This is a bounded invertible operator, 
therefore the inverse, that we denote by Q.~ g , is a bounded operator 
from Q q to Qq. This inverse operator is a compact perturbation of 0. q . 

We will use Proposition 3.4 and Lemma 3.5 many times throughout 
the paper. An immediate useful corollary of the latter lemma is the 
following. 
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Corollary 3.8. Under the conditions of Proposition 3.4, for any e > 
there exists a subspace jj( C Ho such that 

\\Q q v\\ 2 > (l-e^lMI 2 , veJ?, 

and the codimension of M in Ho is finite. 

4. Approximate spectral projections 

The aim of this section is to prove that the subspaces Q q are very good 
approximations to the actual spectral subspaces H q of the operator , 
as well as to the spectral subspaces of ^_(V). Closeness of subspaces 
will be measured by closeness of orthogonal projections onto them. 
Recall that the projection onto H q is denoted by P q . Let Q q be the 
orthogonal projection onto Q q . 

Theorem 4.1. The projections P q and Q q are close: for any N , and 

any collection of the operators Tj, j = 1, . . . , N, each ofTj being or 
0, the operator T(P q — Q q )T , is compact, T = T{T 2 . . . T^. 

In justifying the theorem, we need two technical lemmas which will 
be proved in Section 7. 

Both lemmas concern the operator properties of products of many 
copies of the resolvents of and ty^, the creation and annihilation 
operators and 0, functions hj with all derivatives bounded, and, 
possibly, the spectral projection P q onto H q . In such product, we assign 
order 1 to and 0, order —2 to the resolvent, order to functions and 
projections. The order of the product is defined as the sum of orders 
of factors. 

Lemma 4.2. Let the operator 21 be the product of creation, annihilation 
operators, resolvents, and functions hj, have negative order, and let at 
least one of the functions hj have compact support. Then 21 is compact. 

Lemma 4.3. Let the operator 21 be the product of creation, annihilation 
operators, resolvents, functions and a spectral projection P q . Then 21 
is bounded. If, additionally, at least one of the functions has compact 
support then 21 is compact. 

The proof of the theorem goes in the following way. We construct an 
intermediate operator & q with range in Q q and prove that & q is close 
both to P q and Q q . 

For q > we define the operator & q by the relation 

& q = C- 1 0"P o 9 , C q = q\(2B°y. (4.1) 

So, our expression for the approximate spectral projection is just a 
natural modification of the exact formula (2.7) for the unperturbed 
operator. Equivalently, the operator & q can be described by the for- 
mula 

& q = Cr 1 ^*, (5 = Q q P ; 
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in Proposition 3.4 this operator is shown to be bounded. 

The proof of Theorem 4.1 will consist of two parts, showing that & q 
is close to P q and showing that it is close to Q q . 

Proposition 4.4. The operator & q is close to the projection P q . 

Proof. We give the detailed proof for q — 1 first, and and then explain 
the natural changes one should make for larger q. 

We denote by R-(() the resolvent of the operator we also need 
the resolvent -R+(C) of the operator The projection P can be 

expressed via Riesz integral 

P = (2m)~ 1 [ itL(CR, 

where r is the closed curve, encircling the zero point A , defined, 
together with curves T q , in Sect. 3. We are going to transform the 
expression for the resolvent -R-(C) using the commutation relations 
(2.9), (2.10), (2.11), and (3.8). After this, the crucial observation is 
that the integral of k > 2 along a closed contour vanishes. This 
enables us to dispose of terms that are not compact before integration. 
We start by writing 

itL(o = pp_ - cr 1 = op+ - 2B° - 2b - c)- 1 . 

Next we use the resolvent equation to get 

R_(() = (y + -2B°-()- 1 - (4.2) 

(<p_ - C) _1 (2b)(q3 + - 2B° - C)- 1 = i?+(2B° + C) - # (C) 

We multiply (4.2) by from the left and by from the right, as (4.1) 
requires. For the first term we use that 

or+(2b° + on = 0(00 - 2B° - cr x = (4.3) 

00(013 - 2B° - Cr 1 = «P_i2_(2B° + C). 
Integration gives 

/ 0(^_-2B°-C)- 1 0rfC = ( P- / (qJ_-C)" 1 dC = 27ri*P_P 1 . 

By Theorem 2.1, the point Ai = 2B° is an isolated point of the essen- 
tial spectrum of ^P_, therefore the nonzero eigenvalues of may 
converge only to A b and thus (<£_ - 2B°)Pi = <p_Pi - 2B°Pi is com- 
pact. Moreover, by Proposition 3.3, this compactness is preserved after 
multiplication by any product of creation and annihilation operators. 

Now we consider the second term, H(() in (4.2). We are going to 
show that for any operator of the type T the integral of TH(()T along 
To is compact. First we apply the resolvent formula (4.2) k times to the 
first factor in H((), with k large enough. This operation will produce 
terms of order —4, —6, . . . , containing factors R-(() and b, and the re- 
mainder term of order —2k containing these factors with, additionally, 
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one factor i? + (2B° + (). This last remainder term, for k sufficiently 
large, is compact by Lemma 4.2, even after the multiplication by T 
from the left and from the right. 

The leading terms in H((), having orders —4, —6,. . . , will be trans- 
formed by repeatedly commuting b and R-(C) and then the resulting 
commutants again with -R-(C) and so on. Under commuting -R-(C) 
with a function, with 0, or with 0, this factor R-(() moves to the left 
or to the right, and one more product in the sum composing TH(()T 
arises, of order lower by 1. In this commuting procedure we aim for 
collecting the factors -R-(C) together. As soon as a term with all R-(() 
standing together arises, we leave it alone and do not transform any 
more. After sufficiently many commutations, we arrive at a collection of 
terms in TH(()T of negative order, which are compact by Lemma 4.2. 
The terms of nonnegative order, the ones to which this Lemma can- 
not be applied, will have the form G 1 R_(() k G 2 with k > 1 and some 
operators Gi,G 2 - These terms vanish after integration over r . 

Now we demonstrate in detail how this procedure works. Let T = 1. 
For the case q — 1, which we are considering at the moment, it suffices 
to apply (4.2) twice and write 



We multiply (4.4) by and 0, as (4.1) requires, and then the last term 
takes the form S J R_(C)bi?_(C)b J R + (2B° + £);Q. This is an operator of 
order —2, containing as a factor the compactly supported function b, 
so it is compact for any ( G r\ by Lemma 4.2. As for the first term 
in (4.4), we handle it in the following way (we suppress ( as long as it 
does not cause a misunderstanding) 



(we used here that [R-,L] = — L]R_ for any operator L.) 

Again we multiply (4.5) by and 0. In the first and the second 
terms on the right-hand side all entries of the resolvent are collected 
together, so these terms vanish after the integration; the last one has 
order —4 and is compact by Lemma 4.2. If the operator T is present, 
we continue expanding (4.4) sufficiently long: 



H(() = 2i?_(C)bi? + (2B° + C) = 2itL(C)bitL(C) 
-4 J R_(C)bi?-(C)b J R + (2B° + C). 



(4.4) 



R bR = bR 2 __ - b]R 2 _ = 

bR 2 _ - [<p_, b]R 3 _ + R-[¥-, [?p_, bp 3 . 



(4.5) 



H(Q = 2R_(()bR + (2B° + C) = 2R_bR 
-4i?_bi?_bi? + (2B° + C)... 



(4.6) 



k 



J2(-2R-b) j R- + (-2i?_b) fe+1 i? + (2B° + C). 
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We choose k large enough, so that the last term in (4.6) will be of neg- 
ative order even after the multiplication by T. In all remaining terms 
in (4.6) we commute R_ with other operators, moving all factors i?_ 
say, to the right. For example, the term i?_bi?_bi?_ will be handled 
in the following way 

R^hR_hR_ = hR 2 _ bi?_ + R_\b, yS_]R 2 _bR- 
= bi?_bi? 2 +bi^[b,<p_] J R_bi?_ + J R_[b,<p_]i? 2 b J R_. 

The second and third terms in the last line are operators of order —7, 
and they stay compact after multiplication by four operators 0, 0, 
which gives compactness of corresponding terms in T0if (z)0T for an 
operator T containing just one factor. The first term in the last line 
has order —6, so we have to make one more commutation: 

bi?_bi? 2 = b 2 R 3 _ + bR [b, <$-]R 3 _. (4.7) 

Now the first term in (4.7) vanishes after integration, and the second 
one has order —7. If the operator T contains more terms, we have to 
continue commuting, moving the factors i?_ to the right. If there is a 
separately standing R_ in some term, we can move it to the right, and 
the order of the term with commutant decreases by one. We cannot 
do this and, correspondingly, decrease the order, only if all R_ stand 
together, but such terms vanish after the integration. This reasoning 
establishes our property for & 1 . 

Now we explain how one handles the operators & q , q > 1. It is 
convenient to use the induction in q. So, supposing that for all q < p 
the compactness is proved, we start transforming the expression (4.1) 
for q = p. First we extract the principal term. We write, using (4.2), 

O p R-(0O p = a p (R+(2B° + ()- H(())£l p . (4.8) 

The term H(Q produces a compact contribution, in the same way as 
for p = 1. For the leading term in (4.8) we use (4.3): 

Q P R + {2B° + C)0 P = 5 p " 1< P_i?_(2B° + C)H P_1 - 

We commute £P -1 and in the last expression to get 

~5 p R+ (2B° + ()£T = < p_5 J, - 1 i?_(2B° + C)0 P_1 + (4.9) 

[0 p - 1 , ( P_] J R_(2B o + C)0^ 1 . 

The first term in (4.9) has the form, similar to the one we have started 
with in (4.8), just with p replaced by p—1, and the inductive assumption 
can be applied. The commutator [0 9 ~\ in the second term in (4.9) 
is the sum of p — 1 expressions of the form 

J '[0, ^_]5 q - 2 - J = 2B°0 p - 1 + 20 J + 1 b0 p - 2 - J , j = 0, . . . ,p - 2. 

(4.10) 
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This gives additionally p — 1 terms of the form 2B°£P 1 i?_(2B° + 
C)£P -1 in our integrand. So we have transformed 0. P R+(2B° + C)5 P 
to (2(p - 1)B° + < p_)0 p - 1 i?_(2B° + C)£F _1 plus remainder terms. 
After p such steps we arrive at the expression 

v 

3Ki2_(2pB° + C); OK = J](2(p - j)B° + ?p_), 

which, after the change of variables £' = 2pT5° + ( and integration along 
T p , gives 27ri9JLP p . The latter operator is a compact perturbation of 
2mC p Pp, by Proposition 3.1, and, as usual, this perturbation remains 
compact after multiplication by T. 

The remainder terms arising in the expressions similar to (4.8), are 
taken care of by means of Lemma 4.2, as in the case of q — 1: they 
contain at least two resolvents separated by a function and the com- 
mutation procedure succeeds. Another kind of remainders are the ones 
obtained in the relations similar to (4.10), substituted into (4.9). They 
contain only one resolvent. However after the integration such terms 
take the form T'PT" where P is one of projections P q , and T', T" are 
products of creation, annihilation operators and functions, moreover at 
least one of these functions has compact support. By Lemma 4.3, such 
operators are compact. □ 

So we have proved that the operator & q is close to the projection P q . 

Proposition 4.5. The operator & q is close to the projection Q q onto 
the approximate spectral subspace Q q . 

Proof. First of all, since the operator & q is a self-adjoint operator close 
to some projection, it must have the essential spectrum as any projec- 
tion has, this means, consisting only of the points and 1. Moreover, 
since the range of & q lies in Q q , the operator & q has Q q as an in- 
variant subspace and it annuls the orthogonal complement to Q q . So 
we may consider & q in Q q . We are going to show first that zero is 
not a limit point of nonzero eigenvalues of & q ; this will imply that 
zero is an isolated point of the spectrum of & q and the range of & q 
is closed in Q q . To show this, we consider the self-adjoint operator 
C q l <5*<5 : Ho -> Ho, (5 = 9 P O - By Proposition 3.4 and Lemma 3.5, 
this operator is bounded, and, moreover, 

C q l <&*<& = P (1 + £)P (4.11) 

with a compact operator This implies that the point is an isolated 
point of the spectrum of C' 1 ®*® : Tio — > Tio, and therefore is an 
isolated point of the spectrum of & q = C^ 1 (25(25*. So the range of & q is 
closed. Further on, it follows from (4.11) that (25* = P £ 9 : Q q -> 
has range of finite codimension in Ho, and since £l q maps H onto 
Q q (by the definition of the latter space), the range of & q has finite 
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codimension in Q q . These properties, taken together, imply that & q 
is a compact perturbation of the projection Q q acting in the range 
of Q q , & q = Q q (l + &i)Q q , with a compact operator R\. Since the 
range of & q has finite co-dimension in Q q , this relation can also be 
written as & q — Q q = & q &2& q + $ with a compact operator & 2 and 5 
having finite rank. By Lemma 4.3, the operator & q stays bounded after 
multiplication by any collection of creation and annihilation operators, 
and this establishes closeness of & q and Q q . □ 

Propositions 4.4 and 4.5, taken together, prove Theorem 4.1. 

Now we add a perturbation by a compactly supported smooth elec- 
tric potential V(x). We are going to show that under such perturba- 
tion the spectral subspaces 'almost' do not change. This observation 
justifies the usual physical intuition hinting that in the first order of 
perturbation theory only the eigenvalues change, while the spectral 
subspaces change only in the second order. Note that by our construc- 
tion of the approximate spectral subspaces we have shown that this 
physical intuition is wrong in the case of a magnetic perturbation. 

Note, first of all, that the perturbation of by V is relatively 
compact, therefore, again, the spectrum of the operator = ^P_ + 

V between Landau levels is discrete. We can change the contours T q a 
little, so that they do not pass through the eigenvalues of ty-(V). We 
denote by the spectral subspaces of *$-(V) corresponding to the 
spectrum inside T q , by PY the corresponding spectral projections, and 
by BZ(C) the resolvent of <#-(V). 

Proposition 4.6. The projections P q and Q q are close in the sense 
used in Theorem J^.l. 

Proof. We will prove that the projection is close to P q , then the 
result will follow from Theorem 4.1. 

We use the representation of projections by means of resolvents, 

pv = (2th)- 1 / R v _(()d(. 

The similar expression for P q and the resolvent formula give 

P q v -P q = -(2m)- 1 [ R4()VBZ(()d( = (4.12) 

2N-1 

(^r 1 E / (-R-(ov) k R-(()d( 

k=l 

+(2ni)- 1 [ (-R„(()V) N R v _(()(-VR4()) N d(. 

The last term in (4.12) is compact and stays compact after the mul- 
tiplication by the creation and annihilation operators, as soon as N 
is large enough, by Lemma 4.2. With the leading terms in (4.12), we 
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can perform the same procedure as when proving Proposition 4.4. We 
commute the resolvent with V and with the terms arising by commu- 
tation and so on, aiming to collect the resolvents together all the time. 
As a results we arrive at a number of terms of sufficiently negative or- 
der, thus compact before the integration, and terms with all resolvents 
collected together, and thus vanishing after the integration. □ 

5. Spectrum of Toeplitz-type operators 

We move on to the study of the splitting of Landau levels of our 
operators. Similar to [13], [11], [5], the properties of the eigenvalues of 
the perturbed operators depend on the properties of the spectrum of 
certain Toeplitz-type operators. 

Usually, by Toeplitz operator one understands an operator of the 
form % w = PWP where P is the orthogonal projection onto some 
subspace Q in L 2 and W is the operator of multiplication by some 
function. Alternatively, if we consider the Toeplitz operator in Q the 
operator can be written as % w = PW. 

In [13], [11] such operators were considered, with Q being one of 
Landau subspaces and the electric potential V acting as W. Later, 
in [5], Toeplitz type operators at the lowest LL were considered, with 
some differential operators acting as W. We will study this latter kind 
of Toeplitz type operators, with the space 7io of zero modes of the 
perturbed Pauli operator acting as Q. 

We recall that the projection onto the space TCq = Go of zero modes 
of is denoted by P . We consider linear differential operators £ of 
order / with bounded coefficients, £u = Yla+pKiPapd^d^u. We suppose 
that all coefficients p a p have compact support in the disk \x\ < R and 
I Pa/3 1 < M. With any such operator C we associate the Toeplitz-type 
operator in TLq\ 

% z v = P £*£v, v G H . 
This operator corresponds to the quadratic form 

t £ [ M ] = J \£u\ 2 dx, u e H . 

We are going to study the spectrum of X £ . We denote by = 
A^(£) positive, resp., negative eigenvalues of the operator r £ & . The 
distribution functions n±(\), A > 0, are defined as n±(\) = n±(\] £) = 
#{n : ±A^(£) > A}. The singular numbers s n (£) of the operator 
are just the absolute values of A^ ordered non- increasingly, and their 
distribution function equals n(A; £) = n + (A; £) + n_(A, £). We denote 
by 5(A) the function 

Proposition 5.1. For an operator £ with bounded compactly supported 
coefficients, 

limsupn±(A;£)S(A)" 1 < 1. (5.1) 
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The estimate (5.1) means that the eigenvalues of T £ converge to zero 
extremely rapidly, super-exponentially: 

limsup(n!A^(£))™ < oo. (5.2) 

The estimates (5.1), (5.2) do not exclude the possibility that the opera- 
tor T £ has only a finite number of positive and/or negative eigenvalues, 
counting multiplicity. 

Proof. We start by proving a special version of the proposition, deal- 
ing with the case when the operator £ is just the multiplication by a 
bounded compactly supported function W . The corresponding Toeplitz 
operator will be denoted by % w 

We apply the classical variational approach. Consider the quadratic 
form t w [v] of the operator % w : 

t w [ v ] = J W(x)\v(x)\ 2 dx, v e H . 

The variational description of eigenvalues (Glazman lemma) implies 

n + (X;W) = min codim {J^C H : i w [v] < X\\v\\ 2 , v G JSf}. 

We recall the description of the space Tio to get 

n + (X;W) = min codim pTc F h : t W {e~*f] < X\\e~^ f\\ 2 , f E Jt), 

(5.3) 

where jF b is the space of entire analytical functions, defined in Sect. 2. 

Now we are going to estimate the quantity in (5.3) from above by 
means of changing the forms and the spaces. First, if we replace 
II e ^f II 2 i n (5-3) by something smaller, then there will be fewer sub- 
spaces ^# where the inequality holds and therefore min codim may be- 
come only larger. So we can replace He - */!! 2 by C||e~( 1+e ' >,1 ' /|| 2 where 
\1/ is the potential defined in (2.3) and C is some constant, since, due to 
our estimates of discussed in Sect. 3, e~* > Ce _ ' 1+ ^ $ °. After this, 
if we replace the enveloping space jF b in (5.3) by a larger one, jF( 1+e ) B o, 
then the codimension of subspaces may only grow. Therefore, we arrive 
at 

n+(\,W)< (5.4) 

min codim {JC C : t w [e~*f] < A||Ce _ ^ 1+ ^*°/|| 2 , / e JC). 

Now we notice that t w [e"*/] = J We-™\f\ 2 dx, and We~ 2 * is again 
a bounded function with compact support. Therefore we can apply 
estimates for the right-hand side in (5.4) established in [13], Proposi- 
tion 3.2 or [11], Lemma 6.1, which lead to the inequality (5.1). 

The estimate we have just established can even be improved, if we 
use the results of [5] to majorize the right-hand side in (5.4): 

i B° 

limsup(n!A^(Vy))« < — Cap (supp W), 
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where Cap (supp W) is the logarithmical capacity of the support of W. 

Now we pass to the general case. Let all coefficients p af 3 of the 
operator £ have support in the disk with radius R centered at and 
\p*fi\ < M. 

Denote by \ the characteristic function of the disk D of radius 2R. 
By Cauchy integral formula, applied as in Lemma 3.5, each term of the 
form papdfd^u, u E Ho can be expressed as T a p(x u )t with a bounded 
operator T a p. The upper bound for the norm of this operator depends 
only on the order of derivatives a + (3, on the size of the coefficient p^ 
and on the radius R. Since the quadratic form t £ of the operator T £ 
has the form t £ [u] = || J2Pai3dfdiu\\ 2 = || Y^Papdf <9fx«|| 2 , we have 

|t £ M < x^WxuW 2 = x _1 t x [w], ueH (5.5) 

for a certain x. Thus the singular numbers of the operator are 
majorated by the ones of the operator % w with W = x-, and for such 
operators the required estimate is already proved. □ 

It is more convenient for our needs to formulate Proposition 5.1, 
directly in the variational form, declaring the existence of subspaces 
with prescribed properties. 

Proposition 5.2. There exists a constant x depending only on I, R, M ; 
such that for any e > and for any A > 0, sufficiently small, there ex- 
ists a subspace G Ti° such that the codimension of ^#(A) is no 
greater than (1 + e)H(A) ; and for any order I differential operator £ 
with coefficients p a p supported in the disk with radius R and satisfying 
\p*p\ < M 

t £ [w] < xA||w|| 2 ,w G Jt{X). (5.6) 

The important feature in this formulation is that the subspace ^#(A) 
of controlled codimension can be chosen in such way that it services 
simultaneously all Toeplitz-like operators generated by differential op- 
erators £ of fixed order and with coefficients subject to the above re- 
strictions. 

Proof. By (5.5), as soon as the chosen subspace ^(A) services the 
Toeplitz operator T x , the inequality (5.6) holds for all operators £ 
subject to our conditions. □ 

We formulate the eigenvalue estimate from below for operators 1 W 
only, under the condition that the function W is non-negative. This 
latter condition can be somewhat relaxed, however some kind of pos- 
itivity requirements still remain. In fact, starting from [13], [11], the 
positivity condition for lower spectral estimates remains a serious ob- 
stacle in this field. 
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Proposition 5.3. Let the function W be non-negative and be greater 
than some c > on an open set. Then 

liminfn + (A;W^)S(A)" 1 > 1. 

Proof. We use the variational description of n±(X; W): 

n+(X;W) = maxdim{i? C H : toH < A|M| 2 , v E &,v ^ 0}. (5.7) 

Again, using the description of the space Ho of zero modes, we can 
rewrite (5.7) as 

n+ (\; W) = maxdim{^ C F h : to[e - */] > Ml^ff, f e ^, f ^ 0} 

(5.8) 

By replacing ||e */|| 2 by some larger quantity, we narrow the choice 
of subspaces ^ where the inequality in (5.8) holds. So, we replace 
lie"*/!! 2 by WCe^-^fW 2 with an arbitrary e 6 (0,1) and some C. 
This leads to decreasing of the right-hand side of (5.8). A further 
decreasing is obtained by narrowing the enveloping space: we replace 
jF b by -7-(i_e)B°- This transforms the right-hand side in (5.8) to the 
form, where the estimate from below is already proved, again in [13], 
Proposition 3.2 or [11], Lemma 6.1. □ 

We will need the uniform variational version of the latter estimate. 

Proposition 5.4. Let Q be an open set, c > 0. Then for any e > and 

for any A > ; small enough, depending on e there exists a subspace 
& = -S?(A) C Ho such that (Wv,v) > \\\v\\ 2 , v e & \ {0} for all 
functions W > satisfying W > c on Q and dimJzf(A) > (1 — e)S(A). 

The result follows from the estimate (Wv,v) > c(xo.v-,v). 

6. Perturbed eigenvalues 

The passage from the spectral estimates for Toeplitz-type operators 
to the ones for the Pauli operator has been performed in a different way 
in [12], [13] and in [11]. A short and neat reasoning in [12], [13] based 
upon the proper version of the Birman-Schwinger principle requires 
positivity of the perturbing operator which fails to be positive in our 
case. Therefore we follow the approach of [11], Sect. 9, using variational 
considerations, however certain new features appear. 

The particular form of the variational principle we use, derived im- 
mediately from the spectral theorem, is the following. For a self-adjoint 
operator 21 and an interval (s, r) in the real axis we denote by N(s, r; 21) 
the total multiplicity of the spectrum of 21 in (s, r); when this quan- 
tity is finite, it is the number of eigenvalues of 21 in (s,r), counting 
multiplicities. 
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Proposition 6.1. Let 21 be a self-adjoint operator with domain 2) (21). 
Then, for s < r , 

N(s,r\%) 

= maxdimj^C £(2t) : || (21 - /i)m|| 2 < T 2 \\u\\ 2 ,ue £\{0}} (6.1) 
= mincodim{j^cD(2l) : \\(%-fi)u\\ 2 >T 2 \\u\\ 2 ,ue^} , (6.2) 

where /i — (r + s)/2 and r = (r — s)/2. 

When using Proposition 6.1, one looks for a subspace where the 
inequality in (6.1) is fulfilled. The dimension of this subspace gives an 
estimate from below for N(s,r; 21). On the other hand, having found 
a subspace where the inequality in (6.2) is satisfied, we can be sure 
that the codimension of this subspace estimates N(s, r; 21) from above. 

Considering perturbations of the Landau Hamiltonian by general 
electric and magnetic fields, we should accept the possibility that the 
eigenvalues may split away from the Landau level both up and down. 
Therefore we need rather advanced notations. We fix a Landau level 
A = A q = 2qB° and set s± = A±B°. For given A G (0, B°/4) we define 
fj,± = (A ± A + s±)/2, r = |B° — |. Thus, /x± is the midpoint of the 
interval between A ± A and s±, r is the half-length of this interval. For 
the perturbation given by an operator 2J we denote by N + (\; q, b, QJ) 
the eigenvalue counting function N(A q + \, s + ; *P_ + 2J) where is the 
Pauli operator with magnetic field B = B°+b; similarly, N_(X; q, b, QJ) 
is defined as iV(s_, A q — A; ^3_+5J). Thus, A^ + , AL count the eigenvalues 
converging to A from above, resp., from below. 

6.1. Upper estimate. We start with the upper estimate. It declares 
that the eigenvalues split away from the Landau level at least superex- 
ponentially rapidly. This estimate is rather robust: is does not require 
much smoothness of the perturbing magnetic field or of the electric 
potential. It is convenient to keep all the previous notations in the 
paper for a smooth case, and use the superscript ft for the non-smooth 
versions. 

So, let h$ e Cq(M 2 ) be a magnetic field with compact support, ^ G 
C 4 be its scalar potential, A^" = b", and a 11 e C 3 its vector potential. 
The Pauli operator associated to the field B° + b" will be denoted by 

*pl. 

We consider the smoothened magnetic field, b = u * h^, where 
uj G is a smooth function with compact support, f udx = 1, and 
u(x) depends only on \x\. Thus the field b is smooth and compactly 
supported. Of course, ip = u * ip^ is a potential for b, Aip = b. More- 
over, since is a harmonic function outside the support of b", by 
the mean value property, ip coincides with ^ outside some compact. 
Correspondingly the vector potential a of the field b coincides with a" 
outside a compact. 
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As a result, the Pauli operator 53_ can be expressed as 

q£ = <p_ + 3JI, 

where 53 _ is the Pauli operator with smooth magnetic field B = B° + b 
and 3JT is a first order differential operator with C 2 -coefncients having 
compact support. 

Further on, we consider the operator 5^1 + V as a perturbation of 

+ v = <p_ + 53, 53 = an + v; v e c 2 (r 2 ). 

Now we can formulate the main upper estimate. 
Theorem 6.2. If the perturbations b" and V belong to Cq(R 2 ) then 
limsupA^ ± (A,^L + \/)S(A)- 1 < 1. (6.3) 

Proof. We will prove the 'minus' version of the Theorem; the 'plus' 
version differs just by non-essential details. Keeping in mind Proposi- 
tion 6.1, we are going to construct a subspace J*?in the domain of the 
operator 53 _ such that the inequality 

||(<P_ + 9J-^)u|| 2 > r 2 |M| 2 

holds for all u £ Jzf , and estimate its codimension. 

We represent an arbitrary function u as u = u\ + u 2 where u\ £ 
Q = Qq belongs to the approximate spectral subspace constructed in 
Section 3 and u 2 belongs to the subspace Q orthogonal to Q, u x = Q q u, 

U 2 = QgU, Q q = l- Q q . 

Then the inequality we aim for takes the form 

||pp_ + 53 - + u 2 )\\ 2 - r 2 \\ Ul \\ 2 - r 2 \\u 2 \\ 2 > 0. (6.4) 

The left-hand side in (6.4) can be written as 

[||(qj_ + _ ^) Ul f - T *\\ Ul f] + [||(qj_ + - v)u 2 \\ 2 - r 2 \\u 2 \\ 2 } 

+23?((53- + 53 - n) Ul , (53- + 53 - n)u 2 ). (6.5) 

Following the pattern of the proof of Proposition 9.2 in [11], we are 
going to find subspaces of controlled codimension in Q and in Q, where 
the first and the second terms in (6.5) are positive, with some margin, 
while the third term in (6.5) is majorized by this margin. 

The first term. We start with the first term. Since U\ £ Q, by our 
construction, u\ can be written in a unique way as u\ = Q. q vi, v± £ Tio, 
and thus the term takes the form 

||(53_ + 53 - ^K|| 2 - r 2 |M| 2 = ||(<p_ + 53 - ^)5%|| 2 - r 2 ||5%|| 2 

(6.6) 

= ||(53_ -/x)S^i|| 2 -r 2 ||S^i|| 2 + 2^((q3_- A t)S^i,53S^i) 

+(53n%,530%). 
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The first two terms in (6.6) can be transformed by means of Proposi- 
tion 3.4; they produce C q ((A — fi) 2 — r 2 )||i;i|| 2 as well as C q (WiVi,Vi) 
with some compactly supported function W\. The last term is non- 
negative and can be ignored when estimating the whole expression 
from below. As for the remaining, next to last term in (6.6), we have 

2K(C£_ - /i)5%, QJ5%) > ~(\x^vxf + ll^if ), 

where x&i, &2 are order q + 2 and order q + 1 differential operators 
with bounded compactly supported coefficients, £i = (*$_ — fi)£l 9 , 
£2 = and x is a characteristic function of an open set containing 

the support of coefficients of 03, so that x^ — 
As a result, we obtain 

m _ + 5J _ ^) Ul f _ r 2|| Ml ||2 > C?((A _ ^2 _ r 2 )||wi|| 2 (67) 

The expression ((A — jj) 2 — r 2 ) equals (A — ji — r)(A — /x + r) = 
A(A — s). Now we are going to chose our first subspace. We fix e G 
(0, |). By Proposition 5.2, there exists a subspace M\ in 7i such that 
|(WVl, t>i)|, ||x£i^i|| 2 , H^^ill 2 are all no greater than ^jA(A — s)||t>i|| 2 
and the codimension of in Ti.Q is less than (1 + e)S(^(A — s)A) 
for A small enough. So, for v\ G the right-hand side of (6.7) 
is no less than \C q (A — s)A||t>i|| 2 . This means that on the subspace 
jSfi = which has the same codimension in Q as M\ has in 

Ho, the left hand side of (6.7) is no less than \C q (A — s)A||t>i|| 2 . The 
operator 9 : 7i — > is boundedly invertible, so we can majorize 
1 1 fill by I j 1*1 1 1 , as in Corollary 3.8, therefore, for some constant Ci, not 
depending on e, A, 

||(q3_ + 2J-//)wi|| 2 -r 2 ||wi|| > ciAHmxII 2 , Ul G Jzfi. 

The subspace jSfi has codimension no greater than (1 + e)S(^(A — s)A) 
in 

The second term. In the next term in (6.5), ||(^3_ — pi + 2J)m 2 || 2 
-r 2 ||-u 2 || 2 , we can write u 2 = (l-Q q )u 2 = (l—P q )u 2 +(P q —Q q )u 2 where 
P q is the spectral projection of corresponding to the neighborhood 
S q — (A — 7, A + 7) of A = A q , defined in Section 3. We denote by P 
the spectral projection of corresponding to the set (s — B°/2, A — 
7] U [A + 7, A + B°). Since the spectrum of is discrete in (A — 
2B°, A) n (A, A + 2B°), the projection P has finite rank, and thus u 2 
can be expressed as 

u 2 = (1 - P q - P)u 2 + (P q -Q q + P)u 2 = g + h. (6.8) 

Constructed as above, the function g is lying in the spectral subspace 
of corresponding to the exterior of the interval (A — 3B°/2, A + B°). 
So, by means of the spectral theorem, we can estimate ||(^3_ — n)g\\ 
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from below as ||(^P- — fJ-)g\\ > ij + B°/2)||p||. Therefore 

||(^--^ir-r 2 y| 2 >c 2 y| 2 , (6.9) 

as well as 

||(q3--^|| 2 -r 2 y| 2 > C3 ||(^_-^|| 2 (6.10) 

with some constants c 2 ,cs not depending on A. 

Now let us consider h = (P q — Q q + P)u 2 . As we have proved 
in Section 4, the operator P q — Q q is compact and remains compact 
after multiplication by any combination of creation and annihilation 
operators. This property remains intact after adding the finite rank 
operator P. So, we can write (6.8) as 

u 2 = g + &u 2 (6-11) 

with a compact operator Therefore, on a subspace of finite codi- 
mension, (6.11) can be rewritten as 

u 2 = g + &g 

with a compact operator 8! which, again, remains compact after mul- 
tiplication by 0,0. As a result, for e fixed, we can find a subspace of 
finite codimension in Q 1 n Dom(<p_) such that 

IM e(l-e,l + e), ll£--^ll e( i- e> i + e ). 



IMI v ' " \\W--n)9\\ 

Therefore the inequalities of the form (6.9), (6.10) can be written for 
u 2 as well, just with a slight worsening of constants, the value of which 
is not important: 

H0p_ -/i) M2 || 2 -t 2 imi 2 > c 4 ||m 2 || 2 , 

as well as 

||(qj_ - ^ 2 || 2 - r 2 |M 2 > c 6 ||(«p_ -fi)u 2 \\ 2 . (6.12) 

Now we can consider the complete second term in (6.5). This ex- 
pression can be written as 

HOP- + 2T - »)u 2 \\ 2 - r 2 \\u 2 \\ 2 = ||(qj_ - n)u 2 \\ 2 - r 2 \\u 2 \\ 2 (6.13) 

+ ||Q3m 2 || 2 + 23?((<p_ - fi)u 2 , Wu 2 ), u 2 E <3. 

The form ||5J-u 2 || 2 is compact with respect to ||(^3- — fi)u 2 \\ 2 , therefore 
on a subspace of finite codimension we have ||2Ju 2 || 2 < ^11 OP- - aO^II 2 - 
Further on, by Cauchy inequality, 

|23?(0P_ - n)u 2 ,m 2 )\ < ^||(<P_ - fi)u 2 \\ 2 + -\\Xu 2 \\ 2 . (6.14) 

4 c 5 

Again, using the same relative compactness, we can restrict ourselves 
to a subspace of finite codimension so that the last term in (6.14) is no 
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greater than — ji)u 2 \\ 2 . After substituting these estimates into 

(6.13) and recalling (6.12), we obtain 

||0P_ + - fi)u 2 \\ 2 - r 2 \\u 2 \\ 2 > - v)u 2 \\ 2 > c 6 \\u 2 \\ 2 (6.15) 

on a subspace ££ 2 of finite codimension in Q' D Dom(*p_). Note that 
the subspaces here do not depend on the value of A. 

The third term. We are going to show now that the third term in 
(6.5) is majorized by the sum of two first ones. We can write it in the 
form 

29fc((<p_+9J-Ai)ui, (q}_+9J-Ai)u 2 ) =28fc((<p_+9J-Ai) 2 ui,U2)- (6.16) 

The function u\ G Q q can be written as U\ = 9 t>i for some v\ G Ho- 
So we have 

(<P_ + 23 - fi) 2 Ul = (00 - pif&V! + (<P_ - a/)9Jui+ (6.17) 

23(<£_ - fi) Ul + = (00 - ^) 2 0% + G. 

We start with the first term in (6.17). 

In the expression 000 9 t> 1 we perform our usual commuting pro- 
cedure moving to the right all the time. As soon as reaches the 
utmost right position, the corresponding term vanishes. Thus we arrive 
at the representation 

000^ = 2 ^0 j (B° + b)5 9_ V (6.18) 
i 

= 2qB°Q q v 1 + ^2'O j hQ q - ] v = Am + Fi. 

In a similar way, we transform the expression 0000-ui = 00000 <? t> :L : 
we move both copies of the operator to the utmost right position, 
where they vanish, being applied to V\ G TLq. What is left has the form 

00000% = AV + F 2 . (6.19) 

The functions Fi,F 2 in (6.18), (6.19) are compositions of several, no 
more than q, operators and at least one function with compact sup- 
port, applied to v ± . The term G = ((<£_ -/i)2J + 23(<P_ -fi) + Q3 2 )0«t; 1 
in (6.17) involves the differential operator of order q + 2 with bounded 
compactly supported coefficients, applied to v 1 . So we have 

(<p_ + V - /i) V = (A - /i) V - + F 2 + G. 

We substitute this expression into (6.16) and use the fact that u\ and 
u 2 are orthogonal. This leads us to 

2|K(0P_ + V - /i)V, w 2 )| = 2\&((-2 f iF 1 + F 2 + G),u 2 )\ 

< M\\ -2 f iF 1 + F 2 + G\\ 2 + M- 1 \\u 2 \\ 2 , 

where M > can be chosen arbitrarily. We fix M larger than 16c^" 1 , 
where C5 is the constant in (6.15). After this, we take e smaller than 
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By Proposition 5.2 we can choose the subspace C Ho, con- 
structed when considering the first term in (6.4), in such way that 
the codimension of in Ho is no greater than (1 + e)S(xA) for 
some k and that || — 2/iFi + F 2 + G|| 2 is no greater than yA||vi|| 2 . 
With all parameters chosen in this way, for U\ G S£\ = O q ^i and 
u 2 G Jtf 2 the third term in (6.5) is majorated by the sum of the first 
two terms, and thus the whole expression in (6.5) is non-negative for 
m = Mi + «2 G Jz?i © -5?2- The subspace Jzf = Jzfi © «if 2 has codimen- 
sion no greater than (1 + e)H(/tA) + N(e). As a result, due to e being 
arbitrary and to the fact that S(«A) = 2(A), A — > 0, we have 

limsupA_(A;g,b,5J)S(A)~ 1 = limsup AL(A; q, b, QJ)S(xA)" 1 < (1+e) 

A^O A^O 

Finally, we let e — > 0, and obtain (6.3). □ 

6.2. Lower estimate. Now we establish an estimate for N±(\) from 
below. The conditions for this estimate to hold are expressed in the 
terms of an effective potential constructed for each LL from the mag- 
netic field b and the electric potential V. The complicated form of this 
effective potential reflects the complicated character of the influence of 
the magnetic perturbation onto the behavior of eigenvalues: we remind 
again that the latter dependence is not monotone. Our construction 
requires considerable smoothness of b and V; in order to avoid extra 
technicalities, we suppose here that both are infinitely smooth. 

For a fixed LL A = A q , with the same s±,/j±,t as before, we define 
the effective potentials as 

W± = W±[b,V,\] =-(A±A + 2B°)(s + 2B°)Z,[b] - Z q+2 [b] (6.20) 

+2(/i± + 3B°)Z q+1 [b] - X q [b, 4(2B - b + /i ± )b - (4B - 2/i± + V)V\ 

-2X q+1 [b, V - 3b] - 490>,[b, dV - 2db], 

where the expressions Z, X, y are described in Proposition 3.4. 

Theorem 6.3. Suppose that for A close to zero the potential W± is non- 
negative and, moreover, there exists an open set Q where W± > c > 
for c independent on A. Then 

liminf N ± (X,q,b,V)E(X)- 1 > 1 (6.21) 

A^O 

Proof. We consider the eigenvalues below A and prove the estimate 
(6.21) with 'minus' sign. The other case is proved in an almost identical 
way. 

We are going to construct the subspace S£ in Q = Q q , having dimen- 
sion asymptotically greater than (1 — e)S(A) such that the inequality 

||(<P_ + V-i2)u\\ 2 < t 2 \\u\\ 2 (6.22) 

is satisfied for all u E Jz? \ {0}. To do this, we set u = £l q v, v G Ho, 
and thus reduce both parts of (6.22) to quadratic forms on the lowest 
Landau level using Proposition 3.4 and then apply Proposition 5.3. 
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So, we have 

|| + V - fi)u\\ 2 - t 2 \\u\\ 2 = (<$_u, <$_u) + 2U(Vu, <P_m)+ (6.23) 
(/i 2 - r 2 ) |H| 2 + (Vu, Vu) - 2{i(Vu, u) - 2/i(<p_w, u). 

We substitute u = 9 t>, v e TCo into (6.23) and consider each term. 
First, by (3.5), we have 

(/i 2 - r 2 )|M| 2 = ^ - r 2 )C q \\v\\ 2 + ~ r 2 )(2 q [b]v, v). (6.24) 

By (3.6) applied for U — V and then for U = V 2 , we obtain 

-2fi(Vu, u) = -2fi(X q [b, V]v, v); (Vu, Vu) = (X q [b, V 2 }v, v). (6.25) 

For the last term in (6.23), we have 

-2/i(<p_M, u) = -2n(QQu, u) = -2/i(00w, u) + 4/iB°(w, u) (6.26) 
+4/i(bu,u) = — 2fiC q+ i\\v \\ 2 
-2/i(Z q+1 [b]v, v) + 4/iB°C g |M| 2 + ^B°(Z q [b]v, v) + Afi(X q \b, b]v, v). 

Further on, commuting 0, 0, we get 

2$l(Vu, qj_u) = 2$t(Vu, 00m) = 2$t(QVu, 0m) - 4(BVw, u) (6.27) 
= 2(VQu,Qu) +23?([0, V>,5u) - A(BVu,u). 

To the first and the third terms in (6.27), we apply (3.6), and to the 
second term we apply (3.7) (with U = —2idV = [0, V]), obtaining 

2$t(Vu,^„u) (6.28) 
= 2(X q+1 [b, V]v, v) - A(X q [b, BV]v, v) + 43(^[b, dV]v, v). 

Finally, for the term (*P_w, <J}_-u) we have 

||^_w|| 2 = ((00 - 2B)u, (00 - 2B)u) = (<#Mu,£lu) (6.29) 
+4B o2 (w,w) +4((2B° + b)bu,u) -4B°||0w|| 2 -4K(0bw,0«). 

All these terms are of the form already considered in (6.26), (6. 27), (6. 28). 
Thus we obtain 

(¥_u,¥_u) = C q+2 \\v\\ 2 + (Z q+2 [b]v,v) -2B°C q+ 4vf (6.30) 
-2B°(Z q+1 [b]v,v) - 2(X q+1 [b,b]v,v) +4B o2 (C (? |M| 2 + (Z q [b}v,v)) 
-4B (C 9+1 |M| 2 + (Z q+1 [b]v,v)) + A(X q [b, (2B° + b)b]v,v) 
-A{X q+1 [b, b]v, v) - 89f(y,[b, db]v, v). 
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Now we collect the expressions (6.24), (6.25), (6.26), (6.28) (6.30) to 
get 

m _ + V -,x)«f -r^uf^ 2 - r 2 ){C q \\v\\ 2 + (Z q [b]v,v)) (6.31) 
-2fi(X q \b, V]v,v) + (X q [b, V 2 ]v,v) - 2fiC q+1 \\v\\ 2 - 2^(Z q+1 [b]v,v) 
+4/iB l7 9 |H| 2 + 4fiB°(Z q [b]v, v) + MX q [b, b]v, v) + 2(X q+1 [b, V]v, v) 
-A(X q [b,bV]v,v) +43f(y q \b,dV]v,v) + C q+2 \\v\\ 2 + (Z q+2 [b]v,v) 
-2B°C 9+1 ||^|| 2 - 2B°(Z q+1 [b]v,v) - 2(X q+1 [b,b]v,v) + 4B o2 l7 9 |M| 2 

+AB o2 (Z q [b}v,v) - AB°C q+1 \\v\\ 2 - AB°(Z q+1 [b]v,v) 
+A(X q [b, (2B° + b)b]v, v) - A(X q+1 [b, B]v, v) - 8%(y q [b, db]v, v). 

The terms with ||v|| 2 have the coefficient 

(/i 2 - r 2 )C q - 2fiC q+1 + 4fiB°C q + C q+2 - 2B°C q+1 + AB o2 C q 
-4B°C, +1 = ((/i 2 - r 2 ) - 2//A + A 2 )C q = A(A - s)C q . 

The sum of all remaining terms in (6.31) equals exactly (— W_v,v), 
where W- is given by (6.20). Now, by Proposition 5.4 it is possible to 
find a subspace Jl{X) C H such that (Wv,v) > \(L — s)C q \\v\\ 2 for 
v E Hq and dimension of ^#(A) is greater than S(A(L — s)C q )(l — o(A)). 
This is the subspace we need. □ 

If the conditions of Theorem 6.3 are satisfied then the conditions 
of Theorem 6.2 are satisfied as well and we arrive at the asymptotic 
formula 

lim N±(\, q, b, V)E(X)- 1 = 1 (6.32) 

A — >0 

7. Proofs of technical lemmas 

Proof of Lemma 3.2. As usual, we may prove the inequality for func- 
tions u in the Schwartz space first, and then extend it to all functions 
for which the right-hand side of (3.2) is finite, by continuity. 

We will prove Lemma by a double induction. For A^ = 1 the state- 
ment is obvious: ||Hw|| 2 = (^3_-u,w), ||Ow|| 2 = (^3+w, u) = (*p_w,-u) + 
2(B-u,w). Now, suppose that for some A^ we have established (3.2) for 
all N < Nq and for all combinations of Tj. The (exterior) induction 
step will consist of proving (3.2) for A" = N . 

To do this, we are going to show first that one can replace by 
on the right-hand side of (3.2) for all A" < N . This will follow as 
soon as we prove the inequality 

(V-U,u) <C(^u,u) + C'(u,u), N<N (7.1) 

(note that the case N = N is included here.) This fact we again prove 
by induction (interior) in N. The statement is obvious for A" = 1. 
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Suppose that we have proved it for all N < JVi < N and we will check 
it for N = N l . In fact, by (2.10), 

= <pf + % (7.2) 

where is a sum of many terms !tH K . Each !lH k is the product of less 
then Nx factors and some factors 2B. In other words, this means 
that 9i K is a product of no more than Aq — 1 operators 0, no more than 
Ni — 1 operators £3 and several functions, with all derivatives bounded. 
Now we transform each term fR K in the following way. We move half of 
the factors 0, in !tH K to the left and the other half of these operators 
to the right, commuting them with functions and using the relations 
(3.8). 

So, many terms will be produced, having the form 

T\ . . . T k FT k+ i . . . T k+m , (7.3) 

where k, m < Aq, each Tj is £3 or £3 and F is a bounded function. Now 
we can estimate each term in (7.3): 

| (Ti . . . T k FT k+1 . . . T k+m u, u) | = | (FT k+1 . . . T k+m u, T fc * . . . T*u) | (7.4) 
<C||FT fc+1 ...r fc+m u|| \\T*...T*u\\. 

The right-hand side in (7.4) can be estimated by the interior and the 
exterior inductive assumptions, since there are less than Ni and less 
than N factors Tj involved. Summing all these estimates over all terms 
generated by we arrive at 

(9tu, u) < C(q3f ^u, u) + C'(u, u), 

which, together with (7.2) leads to (7.1) for iV = N x . Changing places 
of and we obtain the dual inequality 

(#Ju, u) < C7(qj%, u) + C'(u, u), N < N . (7.5) 

Now, the inductive step in the exterior induction goes rather easily. 
For N = N , let us have iV operators Ti, . . . ,T N . Suppose that T N = £3. 
Then we can write | |T\ . . .T N u\\ 2 as | |T\ . . . Tjv_!t>| | 2 with v = £lu. 
By the inductive assumption and (7.1), we have | |T\ . . . T^_ 1 v \ \ 2 < 
C(?p+~ 1 v,v) + C'(v,v), the first term here equals C(Vp N u, u) so we 
arrive at (3.2). In the case of T N = 0, we write the left-hand side 
in (3.2) as 1 17\ . . . T N -±v \ | 2 with v = £lu. By the inductive assump- 
tion, this quantity is majorized by (^3^ _1 f,f) + (v,v) which equals 
(SP+'U, u) + (ty+u, u). Applying (7.5), we again obtain the estimate 
(3.2). ^ □ 

Proof of Lemma 4-2. If 21 satisfies the conditions of the Lemma and 
we commute some entry of the resolvent with a neighboring function 
or the operators £3 or then all resulting operators will again satisfy 
these conditions. So we perform such commutings by the following 
rule. First, move one copy of the resolvent to the utmost right position. 
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Then if there are three or more operators 0, in a row in front of some 
resolvent, commute the utmost right of these operator with resolvent. 
After performing this operation as long as it is needed, we arrive at the 
situation where in each summand no more than two operators 0, 
stand in front of each copy of the resolvent. Each such group is a 
bounded operator, by Lemma 3.2. At least one copy of the resolvent 
has no more than one operator 0, in front of it, because the order of 
the whole product is negative. After this, we move the function h that 
has compact support to the position in front of this particular copy of 
the resolvent. In this process, the structure of all products remains the 
same, moreover, when we commute h with resolvent, a compact term 
arises. Finally, we arrive at the term hTR± or hR±, where T is or 0. 
This operator is, obviously, compact. All other factors in the product 
are bounded, therefore the whole product is compact. □ 

Proof of Lemma 4-3. Boundedness of the operator in question is de- 
rived from Lemma 3.2 in the same way as Proposition 3.3 is proved. 
By means of commuting, preserving the overall structure of all arising 
terms, we arrive at a sum of operators of the form 



where 7^, T K are products of no more than iV operators 0, 0. The 
middle term in (7.6) is bounded by Proposition 3.1, and the extreme 
terms are bounded by Lemma 3.2. If, say, the function h K is compactly 
supported, we can take N larger than the number of terms in T K , and 
then the operator h K T K {ty- + 1)~ N has negative order and is compact 
by Lemma 4.2. □ 
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